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ABSTRACT 

Vega's photospheric abundances of Li, Na, and K were determined by using con- 
siderably weak lines measured on the very high-S/N spectrum, while the non-LTE 
correction and the gravity-darkening correction were adequately taken into account. 
It was confirmed that these alkali elements are mildly underabundant ([Li/H] ~ —0.6, 
[Na/H] ~ —0.3, and [K/H] ~ —0.2) compared to the solar system values, as generally 
seen also in other metals. Since the tendency of Li being more deficient than Na and 
K is qualitatively similar to what is seen in typical interstellar cloud, the process of 
interstellar gas accretion may be related with the abundance anomaly of Vega, as 
suspected in the case of A Boo stars. 
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1 INTRODUCTION 

While Vega (= a Lyr = HR 7001 = HD 172167 = HIP 91262; 
A0V) plays an important role as the fundamental photomet- 
ric standard, its photospheric chemical composition is known 
to be sort of anomalous. Namely, most elements are mildly 
underabundant (by ~ —0.5 dex on the average) in spite of 
its definite population I nature, though the extent of defi- 
ciency is inconspicuous for several volatile elements of low 
condensation temperature (T c ) such as C, N, O, and S. Since 
this is the tendency more manifestly shown by a group of 
A-F main-sequence stars known as "A Boo-type stars" (see, 
e.g., Paunzen 2004 and the references therein), Vega's mod- 
erate abundance peculiarities have occasionally been argued 
in connection with this A Boo phenomenon (e.g., Baschek 
& Slettebak 1998; Stiirenburg 1993; Holweger & Rentzsch- 
Holm 1995; Ilijic et al. 1998). 

Regarding the origin of unusual surface compositions of 
A Boo stars, various models have been propounded so far, 
such as the diffusion/mass-loss model, accretion/diffusion 
model, or binary model (see the references quoted in Paun- 
zen 2004). Among these, the interpretation that the anomaly 
was built-up by the accretion of interstellar gas (where re- 
fractory metals of high T c are depleted because of being 
condensed into dust while volatile species of low T c hardly 
suffer this process), such as the interaction model between 
a star and the diffuse interstellar cloud proposed by Kamp 
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f Based on observations carried out at Okayama Astrophysical 
Observatory (Okayama, Japan). 



& Paunzen (2002), appears to be particularly promising, 
in view of Paunzen et al.'s (2002) recent finding that the 
[Na/H] values of A Boo stars (showing a large diversity) are 
closely correlated with the [Na/H]isM values (sodium abun- 
dance of interstellar matter) in the surrounding environment 
(cf. Fig. 7 therein) . 

Then, according to the suspected connection between 
Vega and A Boo stars, it is natural to examine the photo- 
spheric Na abundance of Vega. Besides, the abundances of 
Li and K may also be worth particular attention in a simi- 
lar analogy, since these three alkali elements play significant 
roles in discussing the physical state of the cool interstellar 
gas through their interstellar absorption (resonance) lines of 
Na i 5889/5896 (Di/D 2 ), Li I 6708, and K I 7665/7699. 

It is somewhat surprising, however, that reliable deter- 
minations of Vega's Na abundance are barely available. To 
our knowledge, although a number of abundance analyses of 
Vega have been published so far over the past half century, 
Na abundance derivations were done only in four studies 
(among which original W\ measurements were tried only two 
of these) according to our literature survey: Hunger (1955; 
Na i 5889/5896), the reanalysis of Hunger's W\ data by 
Strom, Gingerich & Strom (1966), Qiu et al. (2001; Na I 
5896), and the reanalysis of Qiu et al's W\ data by Saffe 
& Levato (2004). Unfortunately, none of these appear to be 
sufficiently credible as viewed from the present knowledge, 
because of neglecting the non-LTE effect for the strong res- 
onance D1/D2 lines which were invoked in all these studies 
(see below in this section). 

The situation is even worse for Li and K, for which any 
determination of their abundances in Vega has never been 
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reported. Gerbaldi, Fraggiana & Castelli (1995) once tried 
to measure the Wa(Li i 6708) value of Vega and derived a 
tentative value of 2 mA, which was not used however, since 
they did not regard it as real under the estimated ambiguity 
of 4 mA (cf. their Table 6). Regarding the K I 7665/7699 
lines, even a trial of measuring their strengths has never 
been undertaken for Vega, as far as we know. 

Several reasons may be enumerated concerning this ev- 
ident scarcity of investigations on Li, Na, and K (which ac- 
tually holds for early- type stars in general): 

— First of all, since the strengths of Li I, Na I, and K I 
lines quickly fade out as T e g becomes higher owing to the 
enhanced ionization of these alkali atoms (only one valence 
electron being weakly bound), one has to contend with the 
difficult task of measuring very weak lines (e.g., W\ value 
even as small as ~ 1mA), which requires a spectrum of con- 
siderably high quality. 

— Second, these lines situate in the longer wavelength re- 
gion of yellow-red to near-IR, where many telluric lines due 
to H2O or O2 generally exist to cause severe blending with 
stellar lines, which often prevent us from reliable measure- 
ments of equivalent widths, (see, e.g., Sect. 6.2 of Qiu et al. 
2001, who remarked the large uncertainty caused by this ef- 
fect in their measurements of Na I D lines). 

— Third, the resonance lines in the longer wavelength re- 
gion are generally known to subject to considerable non- 
LTE corrections especially for the case of Na I 5889/5896 
lines showing considerable strengths. Unfortunately, since 
detailed non-LTE studies of these alkali elements have been 
rarely published as far as early-type (A-F) dwarfs are con- 
cerned □ lack of sufficient knowledge on non-LTE corrections 
may have hampered investigating their abundances. 

Given this situation, we decided to spectroscopically es- 
tablish Vega's photospheric abundances of Li, Na, and K 
scarcely investigated so far, while overcoming the problems 
mentioned above, in order to see how they are compared 
with the compositions of other elements in Vega as well as 
with those of interstellar gas. This is the purpose of this 
study. 

This attempt was originally motivated by our recent 
work (Takeda, Kawanomoto & Ohishi 2007) of publishing 
a digital atlas of Vega's high-resolution (R ~ 100000) and 
high-S/N (~ 1000-3000) spectrum along with that of Regu- 
lus (a rapid rotator to be used for a reference of contaminat- 
ing telluric lines), by which the confronted main obstacles 
(the necessity of measuring very weak lines as well as remov- 
ing telluric blends) may be adequately conquered. Besides, 
regarding the technical side of abundance determinations, 
we are ready to evaluate not only the necessary non-LTE 
corrections for Li, Na, and K (according to our recent ex- 
periences of statistical equilibrium calculations on these el- 
ements for F-G-K stars), but also the gravity-darkening 
corrections which may be appreciable especially for very 
weak lines as revealed from our recent modeling of rapidly- 



1 Few available at present may be only two studies on neutral 
sodium lines: Takeda & Takada-Hidai's (1994) non-LTE analysis 
of various Na I lines in a CMa (in addition to Procyon and A-F 
supergiants), and Andrievsky et al.'s (2002) non-LTE study of 
Na I Di and D2 lines for late-B to early-F type dwarfs of A Boo 
candidates. 



rotating Vega (Takeda, Kawanomoto & Ohishi 2008). There- 
fore, this would make a timely subject to address. 



2 OBSERVATIONAL DATA 

The basic observational data of Vega used for this study 
are the high-S/N (~ 1000-3000) and high resolution (R ~ 
100000) spectra, which are based on the data obtained with 
the High Dispersion Echelle Spectrograph (HIDES) at the 
coude focus of the 188 cm reflector of Okayama Astrophys- 
ical Observatory. These spectral data of Vega, along with 
those of Regulus (a rapid rotator serving as a reference of 
telluric lines) were already published as a digital atlas by 
Takeda et al. (20070 which may be consulted for more de- 
tails. 

The following 11 lines were selected as the target 
lines used for abundance determinations: Na I 5890/5896, 
8183/8195, 6154/6161, 5683/5688, K 1 7665/7699, and Li 1 
6708 (cf. Table 1). Since most of these wavelength regions 
(except for those of Na I 6154/6161 and Li I 6708) are con- 
taminated by telluric lines, these features were removed by 
dividing Vega's spectra by the reference spectra of Regu- 
lus (with the help of IRAF0 task telluric) which are also 
included in Takeda et al. (2007) along with the airmass in- 
formation. Though this removal process did not always work 
very successfully (e.g., for the case of very strong saturated 
lines such as the O2 in the K I 7665/7699 region and the 
H2O lines in the Na I 8183/8195 region), we could manage 
to recover the intrinsic stellar lines more or less satisfacto- 
rily, since none of these target lines fortunately coincided 
with the deep core positions of such very strong telluric fea- 
tures. 

While most of the lines were confidently identified by 
comparing them with the theoretically simulated spectra, 
the detection of the extremely weak Li I 6708 line was still 
uncertain in the original spectrum even at its very high S /N 
ratio of ~ 2000. Therefore, an extra smoothing was applied 
to the spectrum portion in the neighborhood of this line by 
taking a running mean over 9 pixels, by which the feature 
became more distinct (cf. Fig. 6) thanks to the further im- 
proved S/N ratio by a factor of 30 

The measurements of the equivalent widths (with 
respect to the local continuum level specified by eye- 
inspection) were done by the Gaussian-fitting technique, and 
the resulting W\ values are given in Table 1. The relevant 
spectra at each of the wavelength regions are shown in Figs. 

1 (Nai 5890/5896), 2 (Na 1 8183/8195), 3 (Na 1 6154/6161), 
4 (Na 1 5683/5688), 5 (K 1 7665/7699), and 6 (Li 1 6708). 

2 While the spectrum atlas is presented as the elec- 
tronic tables of this paper, the same material is avail- 
able also at the anonymous FTP site of the Astronomical 
Data Center of National Astronomical Observatory of Japan: 
<ftp://dbc.nao.ac.jP/DBC/ADACnew/J/other/PASJ/59.245/) I 

A IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which is operated by the Association of Universities for 
Research in Astronomy, Inc. under cooperative agreement with 
the National Science Foundation, USA. 

4 Though this results in a blurring of ~ 10 km s _1 corresponding 
to the width of 9-pixel boxcar function (1 pixel ~ 0.03 A), this 
does not cause any serious problem because it is smaller than the 
intrinsic line width of ~ 20-30 km s — 1 . 
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Figure 1. Spectra of 5885— 5900 A region comprising Nal 5890 
and 5896 lines. Theoretical simulation (LTE, broadened with 
i; e sini = 22kms _1 ), Vega/Regulus (telluric lines removed), raw 
Vega, and raw Regulus are arranged from top to bottom with ap- 
propriately chosen offsets. The scale of the ordinate corresponds 
to the second-top spectrum, for which equivalent-widths of these 
lines were measured by the Gaussian fitting method (as overplot- 
ted in solid lines on this spectrum). 
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Figure 2. Spectra of 8180-8200 A region comprising Nal 8183 
and 8195 lines. Theoretical simulation, Vega/Regulus (telluric 
lines removed), raw Vega, and raw Regulus are arranged from 
top to bottom with appropriately chosen offsets. Otherwise, the 
same as in Fig. 1. 



3 ABUNDANCE DETERMINATIONS 

Regarding the fiducial atmospheric model of Vega, Kurucz's 
(1993) ATLAS9 model with the parameters of T off = 9550 K 
(effective temperature), log g (cm s -2 ) = 3.95 (surface grav- 
ity), [X/H] = —0.5 (metallicity), and vt = 2 km s _1 (mi- 
croturbulence) was adopted for this study as in Takeda et 
al. (2007), which appears to be the best as far as within 
the framework of classical plane-parallel ATLAS9 models 
according to Castelli & Kurucz (1994). 

As the basic strategy, the LTE abundance (Alte) was 
first derived from the measured W\ by using the WIDTH9 
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Figure 3. Spectra of 6151-6164 A region comprising Nal 6154 
and 6161 lines. Theoretical simulation (upper) and raw Vega 
(lower) are placed with an appropriate offset. This region is not 
contaminated by telluric lines. Otherwise, the same as in Fig. 1. 
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Figure 4. Spectra of 5680-5690 A region comprising Nal 5683 
and 5688 lines. Theoretical simulation, Vega/Regulus (telluric 
lines removed), raw Vega, and raw Regulus are arranged from 
top to bottom with appropriately chosen offsets. Otherwise, the 
same as in Fig. 1. 



program[3, to which the non-LTE correction (Anlte) and 
the gravity-darkening correction (Aqd) were then applied 
to obtain the final abundance (Anlte+gd). 

In this step-by-step wayQ one can quantitatively judge 
the importance/unimportance of each effect by comparing 



5 This is a companion program to the ATLAS9 model atmo- 
sphere program written by R. L. Kurucz (Kurucz 1993), though 
it has been considerably modified by Y. T. (e.g., for treating mul- 
tiple component lines or including non-LTE effects in line forma- 
tion 

6 Complete simulations (where the non-LTE effect and the 
gravity-darkening effect were simultaneously taken into account 
in a rigorous manner) show that the finally derived abundances 
in this semi-classical approach arc fairly reliable. This consistency 
check is separately described in Appendix B. 
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Figure 5. Spectra of 7660-7705 A region comprising K I 7665 and 
7699 lines. Theoretical simulation, Vega/Regulus (telluric lines 
removed), raw Vega, and raw Regulus are arranged from top to 
bottom with appropriately chosen offsets. Otherwise, the same as 
in Fig. 1. 
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Figure 6. Spectra of 6700— 6716 A region comprising Li I 6708 
merged doublet. Theoretical simulation, smoothed Vega spectrum 
by convolving the 9-pixel boxcar function, and raw Vega spec- 
trum are arranged from top to bottom with appropriately chosen 
offsets. This region is hardly contaminated by telluric lines. Oth- 
erwise, the same as in Fig. 1. 



two abundance corrections, and such factorized information 
may be useful in application to other different situations 
(e.g., the case of slow rotator where only Alte is relevant; 
cf. Appendix A). Also, the use of "relative correction" (in- 
stead of directly approaching the absolute abundance solu- 
tion) has a practical merit of circumventing the numerical- 
precision problem involved in the spectrum computation on 
the gravity-darkened model (especially seen for the case of 
the extremely weak Li I 6708 line composed of several com- 
ponents; cf. footnote 9). 

The non-LTE calculations for evaluating Anlte were 
carried out in the same manner as described in Takeda et 
al. (2003; Na i), Takeda & Kawanomoto (2005; Li i), and 
Takeda et al. (2002; K i). Besides, all the relevant atomic 



data (gf values, damping constants, etc.) for abundance de- 
terminations were so adopted as to be consistent with these 
three studies. 

The gravity-darkening corrections (Agd) were derived 
based on the equivalent-width intensification factor com- 
puted by the program CALSPEC under the assumption of 
LTE (W4/W0, where W± and Wo correspond to the gravity- 
darkened rapid rotator model and the classical rigid-rotation 
model, respectively; and the abundance was adjusted to 
make W4 consistent with WjJ ), the details about which 
are explained in Takeda et al. (2008). Practically, the rela- 
tion Agd = — log(W4/Wo) was assumed for the very weak 
lines {W x < 15 mA; i.e., Na I 6154/6161, Na I 5683/5688, 
K 1 7665/7699, and Li I 6708) being guaranteed to locate 
on the linear part of the curve of growth. The Agd val- 
ues for the remaining four lines (Na I 5890/5896 and Na I 
8183/8195) were directly obtained from the abundance dif- 
ference between those derived from the true W^ bs and the 
perturbed W A ° bs /(W4/W0). 

The results of the abundances and the abundance 
corrections are summarized in Table 1, where the finally 
adopted (average) values of (Anlte+gd) and the abun- 
dances relative to the standard solar system compositions 
([X/H] = (^nlte+gd) — A®) ar e also presented. As can 
be seen from this table, the non-LTE corrections are always 
negative (corresponding to the non-LTE line-intensification) 
with appreciably WA-dependent extents (ranging from 0.1- 
0.7 dex). Similarly, the gravity-darkening corrections are 
also negative with extents of 0.1-0.3 dei\j (especially im- 
portant for the weakest Li I 6708 line) reflecting the line- 
strengthening caused by rotation-induced T-lowering. 

Regarding errors in the resulting abundances, those due 
to uncertainties in the atmospheric parameters (e.g., in T e a 
or log g) may not be very significant in the present compar- 
atively well-established case (an abundance change is only 
~ 0.1 dex for a change either in T e ff by ~ 200 K or in 
log 5 by 0.2; cf. Table 5 in Takeda & Takada-Hidai 1994 for 
the case of a CMa). Rather, a more important source of 
error would be the ambiguities in the measurement of Wx- 
While the photometric accuracy of Wx estimated by ap- 
plying Cayrel's (1988) formula (FWHM of w ~ 0.5A, pixel 
size of 8x ~ 0.03A, and photometric error of e ~ 1/2000) 
is only on the order of ~ 0.1 mA and inconsiderable, un- 
certainties in specifying the continuum level can be a more 
serious problem, especially for the case of extremely weak 
lines where the relative depression is only on the order of 



10" 



Presumably, errors in W\ would amount to several 



tens of % in the case of very faint lines (e.g., Na I 6154 or 
Li 1 6708) or the case of insufficient spectrum quality due 



7 It may be worth noting that these moderate Agd values cor- 
respond to Takeda et al.'s (2008) model No. 4 (v e = 175 km s -1 ) 
which they concluded to be the best. If higher-rotating models 
with larger gravity-darkening were to be used, the corrections 
would naturally become more enhanced. For example, if we adopt 
Takeda et al.'s (2008) model No. 8 with v e = 275 km s _1 , which 
is near to the solution suggested from interferometric observa- 
tions (Peterson et al. 2006, Aufdenberg et al. 2006), the extents 
of Agd turn out to increase significantly by a factor of ~ 3 (i.e., 
-0.22, -0.18, -0.31, -0.23, -0.42, -0.43, -0.36, -0.32, -0.34, 
-0.37, and -0.82, for Na 1 5890, 5896, 8183, 8195, 6154, 6161, 
5683, 5688, K 1 7664, 7699, and Li 1 6708, respectively). 
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to incompletely removed telluric lines (e.g., Na I 8183 or 
K I 7665). Accordingly, the abundance results derived from 
such problematic lines may be uncertain by up to ~ 0.2 
dex, which should particularly be kept in mind in discussing 
the final abundances of Li and K being based on only 1 or 
2 lines. 



4 RESULTS AND DISCUSSIONS 
4.1 Na abundance of Vega 

The abundances of sodium could be derived from 8 lines 
of different strengths (2 resonance lines of D1+D2 and 
6 subordinate lines originating from the excited level of 
Xiow = 2.21 eV). The three kinds of abundances (Alte, 
Anlte, Anlte+gd; cf. Table 1) for each of the 8 lines are 
plotted against W\ in Fig. 7. We can read from this fig- 
ure that the serious inconsistency seen in the LTE abun- 
dances is satisfactorily removed by the non-LTE corrections, 
which are significantly WA-dependent from —0.7 dex (for 
the strongest Na I 5890) to —0.1 dex (for the weakest-class 
Na 1 6154/6161/5683/5688 lines of W\ < 10 mA). This may 
suggest that the applied non-LTE corrections are fairly re- 
liable and the adopted microturbulence of Vt = 2 km s _1 
is adequate (which affects the abundances of stronger Na I 
5890/5896 lines; cf. Table 5 of Takeda & Takada-Hidai 1994). 
The final sodium abundance of Vega, obtained by averaging 
Anlte+gd values of 8 lines, is 6.01 with the standard de- 
viation of a = 0.10, indicating a mildly subsolar composi- 
tion ([Na/H] = —0.3) compared to the solar abundance of 
Aq = 6.31. 

As mentioned in Sect. 1, a few previous determinations 
based on D1+D2 lines are available for the Na abundance 
of Vega. Hunger (1955) tentatively derived A ~ 6.7 (coarse 
analysis) and A ~ 7.8 (fine analysis) from VTa (5890/5896) 
= 195/115 mA, though commenting that these abundances 
are very uncertain. Thereafter, based on Hunger's W\ data, 
Strom et al.'s (1966) model atmosphere analysis concluded 
A ~ 7.7(±0.4) (T cff = 9500 K). Then, after a long blank, 
Qiu et al. (2001) obtained A = 6.45 from their measured 
H / a(5896) value of 94 mA, though they also remarked that 
this result is quite uncertain because of blending with wa- 
ter vapor lines. Soon after, Saffe & Levato (2004) derived 
A = 6.37 using Qiu et al.'s W\ data. Accordingly, all these 
previous work reported supersolar or near-solar Na abun- 
dances for Vega, contradicting the conclusion of this study. 
It is evident that they failed to obtain the correct abundance 
from the Na I 5890/5896 lines because of neglecting the non- 
LTE effect, which is substantially important for these strong 
resonance lines. 



4.2 Composition characteristics of Li, Na, and K 

According to Table 1, while potassium is only mildly subso- 
lar ([K/H] ~ —0.2) to an extent similar to sodium (~ —0.3), 
lithium is more manifestly underabundant as [Li/H] ~ —0.6. 
Therefore, we may conclude that the deficiency in Li is 
markedly different from that of Na and K, even consider- 
ing the possible ambiguities of ~ 0.2 dex (cf. Sect. 3). How 
should we interpret these results? Do they have something 
to do with the element depletion (due to dust condensation) 
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Figure 7. Na abundances derived from each of the 8 lines plot- 
ted against the equivalent widths, based on the results in Table 
1. Squares, triangles, and circles correspond to the LTE abun- 
dances, the non-LTE abundances, and the final abundances (in- 
cluding both the non-LTE and gravity-darkening corrections) , re- 
spectively. The horizontal dotted line indicates the average of the 
final abundances (6.01). 



in interstellar gas and its accretion? Since the quantitative 
amount of such a depletion is different from case to case de- 
pending on the physical condition (e.g., see the variety of 
[Na/H]isM in Fig. 7 of Paunzen et al. 2002), it is not much 
meaningful to discuss the "absolute" values of [X/H] here. 
Rather, we should pay attention to their "relative" behav- 
iors with each other or in comparison to those of various 
other elements. 

The [X/H] vcga values (those for Li, Na, and K are from 
this study, while those for other elements were taken from 
various literature) are plotted against the condensation tem- 
perature (T c ) in Fig. 8 (filled circles), where the [X/H]ism 
results of typical interstellar gas in the direction of £ Oph 
(taken from Table 5 of Savage & Sembach 1996; see also 
Fig. 4 therein) are also shown by open circles. The following 
characteristics can be recognized from this figure: 

— (1) The [X/H]vega values, falling on a rather narrow range 
between ~ — 1 and ~ 0, tend to decrease with T c , a qualita- 
tively similar trend to that seen in [X/H] ism- 

— (2) The run of [X/H]v cg a with T c appears to be slightly 
discontinuous at T c ~ 1000 K (i.e., [X/H] Vcg a ~ -0.2 to 
-0.3 at T c < 1000 K and [X/H] Vcg a ~ -0.5 to -0.6 at 
T C > 1000K). 

- (3) Interestingly, the decreasing [X/H]ism with T c also 
shows a discontinuity around this critical T c (~ 1000 K) in 
the sense that the slope of |d[X/H]isM/dT c | becomes evi- 
dently steeper on the T c > 1000 K side. 

- (4) Specifically, the tendency of [Li/H] Vcg a (—0.6) be- 
ing more deficient than [Na/H]v cg a (—0.3) and [K/H]v cg a 
(—0.2) is qualitatively similar to just what is seen in ISM 
([Li/H] ISM < [Na/H] ISM m [K/H] ISM ). 

All these observational facts may suggest the existence 
of some kind of connection between Vega's photospheric 
abundances and those of interstellar cloud, which naturally 
implies that an accretion/contamination of interstellar gas is 
likely to be responsible (at least partly) for the abundance 
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Figure 8. [X/H] values (relative abundances in comparison with 
the solar-system values of Grevesse &; Noels 1993) of various ele- 
ments plotted against the condensation temperature (T c ). Open 
circles indicate the compositions of interstellar gas in the direction 
of £ Oph (cool diffuse clouds), which were taken from Savage & 
Sembach (1996; cf. their Table 5 and Fig. 4). The [Na/H], [K/H], 
and [Li/H] for Vega derived in this study are shown by larger 
filled circles. Besides, Vega's [X/H] values for other elements are 
plotted by smaller filled circles for a reference, which were taken 
from from various sources: Przybilla & Butler (2001) (for C, N, 
and O; cf. their Table 7); Takada-Hidai & Takeda (1996) (for S); 
Adelman & Gulliver (1990) (for Mg, Ca, Ti, Cr, Mn, Fe, and 
Ni); Qiu et al. (2001) (for Si and V). [Note. Since these litera- 
ture [X/H] values of Vega were derived in a conventional manner 
based on classical atmospheric model, further corrections for the 
gravity-darkening effect (Aqd) are to be expected. According to 
Takeda et al. (2008), however, the extents of (negative) Aqd 
(becoming appreciable only for very weak lines and for specific 
elements/stages) are only ~ 0.2 dex at most (e.g., Adelman & 
Gulliver's [Ca/H] value derived from very weak Ca I lines had 
better be reduced by ~ 0.2 dex) while < 0.1 dex in many typical 
cases, which are thus unlikely to cause any significant change in 
the general pattern of [X/H] shown here.] 



peculiarity of Vega, such as being proposed for explaining 
the A Boo phenomenon. This is the conclusion of this study. 
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APPENDIX A: LI ABUNDANCE IN O PEG 

The enhanced deficiency of Li (compared to Na and K) 
was an important key result in deriving the conclusion of 
this study (cf. Sect. 4.2), since we interpreted it as a man- 
ifestation of ISM compositions. In this connection, it may 
be worth mentioning another remarkable very sharp-lined 
Al IV star, o Peg, for which we could also get information 
of the Li abundance. As the only available measurement 
of the Li I 6708 line for early A-type stars thus far to our 
knowledge, Coupry & Burkhart (1992) derived IT A (6708) = 
1.3 mA for this star, which has atmospheric parameters (T e s 
— 9650 K, log g = 3.6) quite similar to those of Vega except 
for its near-normal metallicity ([Fe/H] ~ +0.1; Burkhart & 
Coupry 1991). 

Now, their W\(o Peg) value twice as large as that of 
Vega (0.7 mA) would raise the Li abundance by +0.3 dex. 
Furthermore, since any flat-bottomed shape has never been 
reported in the spectral lines of o Peg (as we can con- 
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Table 1. Atomic data, equivalent widths, and abundance results. 



Species 


RMT 


A 


Xlow 


log gf 


W X 


^4lte 


Anlte 


A GD 


^NLTE+GD 


[X/H] 






(A) 


(eV) 




(mA) 












Nal 


1 


5889.95 


0.00 


+0.12 


133.3 


6.82 


-0.69 


-0.09 


6.03 




Nal 


1 


5895.92 


0.00 


-0.18 


104.2 


6.58 


-0.49 


-0.08 


6.01 




Nal 


4 


8183.26 


2.10 


+0.22 


19.7 


6.16 


-0.22 


-0.12 


5.82 




Nal 


4 


8194.82 


2.10 


+0.52 


45.1 


6.39 


-0.28 


-0.10 


6.01 




Nal 


5 


6154.23 


2.10 


-1.56 


0.8 


6.36 


-0.10 


-0.12 


6.13 




Nal 


5 


6160.75 


2.10 


-1.26 


1.6 


6.36 


-0.10 


-0.14 


6.12 




Nal 


6 


5682.63 


2.10 


-0.67 


4.3 


6.19 


-0.09 


-0.13 


5.97 




Nal 


6 


5688.21 


2.10 


-0.37 


8.4 


6.20 


-0.09 


-0.12 


5.99 






















6.01 


-0.30 


Ki 


1 


7664.91 


0.00 


+0.13 


13.3 


5.23 


-0.28 


-0.13 


4.82 




Kl 


1 


7698.97 


0.00 


-0.17 


10.7 


5.42 


-0.27 


-0.14 


5.01 






















4.91 


-0.22 



Lii 1 6707.756 0.00 -0.43 0.7 3.15 -0.15 -0.26 2.74 

6707.768 -0.21 

6707.907 -0.93 

6707.908 -1.16 

6707.919 -0.71 

6707.920 -0.93 



2.74 -0.57 



Following the atomic data of spectral lines (species, multiplet No., wavelength, lower excitation po- 
tential, and logarithmic gf value) in columns 1-5, column 6 gives the measured equivalent width. 
The results of the abundance analysis are given in column 7 (LTE abundance; in the usual normal- 
ization of H=12.00), column 8 (non-LTE correction), column 9 (gravity-darkening correction), and 
column 10 (non-LTE as well as gravity-darkening corrected abundance). The finally adopted average 
abundance and the corresponding [X/H] value (= A^ ega — A^) are also given at the bottom of the 
section (in italics; columns 10 and 11), where Grevesse & Noels's (1993) values of 6.31 (Na), 3.31 
(Li), and 5.13 (K) were used as the standard solar-system abundances (A^). Note that the analysis 
of W\ (Li I 6708) was done by synthesizing the six component lines of 7 Li (neglecting the contribution 
of 6 Li; cf. Takeda & Kawanomoto 2005), while all the remaining lines of Na and K were analyzed by 
the single-line treatment as was done by Takeda et al. (2003; Na) and Takeda et al. (2002; K). 



firm from the high-quality spectrum atla^f] of this star pub- 
lished by Gulliver, Adelman, & Friesen 2004), the appli- 
cation of the gravity-darkening correction (~ —0.3 dex for 
Vega) may not be relevant here, which again leads to an 
increase of 0.3 dex. Hence, we have [Li/H] p eg (— [Li/H]vc g a 
+ 0.6) ~ 0.0; an interesting result that the photospheric 
Li abundance of o Peg almost coincides with that of the 
solar-system composition. This result, that [Li/H] (as well 
as [Fe/H]) is deficient/normal in Vega/o Peg, may suggest 
that the mechanism responsible for producing the under- 
abundance of these elements in Vega is irrelevant for o Peg. 
Thus, according to our interpretation, o Peg would not have 
suffered any pollution due to accretion of interstellar gas 
depleted in volatile elements of higher T c . 

Yet, we had better keep in mind another possibility 
that the mechanism causing this Li deficiency in Vega might 
different from that of other metals. For example, this un- 
derabundance could be attributed to some process of en- 
velope mixing (e.g., meridional circulation or shear-induced 
turbulence which are supposed to be more significant as a 
star rotates faster), because Li atoms are burned and de- 
stroyed when they are conveyed into the hot stellar interior 

8 Also available at 

http://www.brandonu.ca/physics/gulliver/ccd_atlases.html . 



(T>2.4xl0 6 K). Namely, since we know that Vega rotates 
rapidly (as fast as v c ~ 200 km s _1 ) while o Peg does not 
so much (at least the gravity-darkening effect is not so sig- 
nificant as in Vega), the underabundance of Li in the atmo- 
sphere of Vega might stem from the rotation-induced mixing 
(which is not expected for slowly rotating o Peg), ff this is 
the case, however, the rough similarity in the extent of de- 
ficiency in Li as well as other metals has to be regarded as 
a mere coincidence, which makes us feel this possibility as 
rather unlikely. 



APPENDIX B: CHECK FOR THE FINAL 
ABUNDANCES: COMPLETE SPECTRUM 
SYNTHESIS 

Our basic strategy for deriving the abundances of Na, K, 
and Li in Sect. 3 was as follows. 

- (1) First, the Alte was derived from W£ bB based on the 
classical plane-parallel model atmosphere. 

- (2) Next, the non-LTE correction (Anlte) was derived 
in the conventional way by using this classical model. 

- (3) Then, the gravity darkening correction (Agd) was 
evaluated from the line-intensification factor Wi TE /Wo TE 
(where the abundance was so adjusted as to satisfy W 7 ^™ ~ 
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W% hs ), which was computed by applying the CALSPEC pro- 
gram (Takeda et al. 2008) with the assumption of LTE to the 
gravity-darkened model 4 and the rigid-rotation model 0. 
— (4) Finally, Anlte+gd was obtained as Alte + Anlte + 
A GD . 

Actually, such a phased approach has a distinct merit 
of clarifying the importance/contribution of two different ef- 
fects (the non-LTE effect and the gravity-darkening effect). 
Besides, from a practical point of view, the necessary amount 
of calculations to arrive at the final abundance solution (such 
that reproducing the observed spectrum) can be consider- 
ably saved. 

However, there is some concern about whether such a 
step-by-step approach (treating the two effects separately) 
really yield sufficiently correct results, because both are ac- 
tually related with each other (e.g., how does the largely 
variable non-LTE corrections over the gravity-darkened stel- 
lar surface differing in T or g play roles? How reliable is 
the gravity-darkening correction derived by neglecting the 
non-LTE effect?). Therefore, it may be worth checking the 
validity of the finally derived abundances (Anlte+gd) by 
carrying out a complete spectrum synthesis including both 
NLTE and GD effects simultaneously. For this purpose, 
the CALSPEC program was modified so as to allow inclu- 
sion of the non-LTE departure coefficients (corresponding 
to the different conditions at each of the points over the 
gravity-darkened stellar surface), and the NLTE+GD pro- 
files (R^ 1 ™) were computed for each of the 9 lines (Na I 
5890, 5896, 8183, 8195, 6f54, 6161, 5683, 5688, K I 7664, 
7699, and Li I 6708) by using the already known Anlte+gd 
values (6.03, 6.01, 5.82, 6.01, 6.13, 6.12, 5.97, 5.99, 4.82, 
5.01, and 2.74, respectively). Such obtained R^ 1 ™ profiles 
are depicted (in thick solid lines) in Fig. 9, where the rele- 
vant three kinds of line profiles [i?Q ,LTE (thick dashed line), 
i?4 TE (thin solid line), and Ro TE (thin dashed line)] are also 
shown for comparison. 

The resulting theoretical equivalent widths (W^ 1 "™) 
computed by integrating R^ LTE (X) are 129.7, 101.3, 18.4, 
43.2, 0.8, 1.5, 4.0, 7.8, 12.8, 10.3, and 0.4 ml for these 9 
lines, respectively. Comparing these values with the observed 
equivalent widths (W^ bs ) given in Table 1, we can see a sat- 
isfactory agreement between W|* LTE and W? bB (typically to 
within several percent except for Li I 670qj), by which we 
may conclude that the approach adopted in Sect. 3 is prac- 
tically validated. 



9 We see a rather large discrepancy for this Li line (Wj* LTE = 
0.4 mA, while W£ hB = 0.7 mA). This must be due to the in- 
evitable numerical errors in 

^NLTE ( L ; y 6708) computed by the 
CALSPEC program (synthesizing the flux spectrum by integrat- 
ing the contributions from a number of finite surface elements; 
each 1° X 1° in latitude and longitude), which becomes apprecia- 
ble especially for this case of Li I 6708 line involving synthesis of 
several component lines. As a matter of fact, we can see a spurious 
wavy pattern in the continuum level outside of this line (see Fig. 
Blk) which indicates the existence of numerical problems (even 
if small). It may be worth stressing that the gravity-darkening 
correction Aqq ~ — log(W4/Wo) is not affected by this numer- 
ical error in W because it is canceled by taking the ratio of W. 
In this sense, our semi-classical approach of establishing the final 
abundance (application of two corrections to the classical LTE 
abundance solution; cf. Sect. 3) is surely advantageous as far as 
this Li I 6708 line is concerned. 
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Figure Bl. Theoretical profiles of the Na, K, and Li lines com- 
puted by the CALSPEC program (developed for synthesizing the 
flux spectrum for a given rotationally-distorted gravity-darkened 
stellar model; cf. Takeda et al. 2008), which has been modified 
to allow inclusion of the non-LTE effect. The ^4nlte+GD values 
given in Table 1 were assumed as the abundances for each of the 
lines. Results for the gravity-darkened model and the classical 
rigid-rotation model (model No. 4 and No. in Takeda et al. 
2008) are discriminated by the line type (solid and dashed lines, 
respectively), and those for NLTE and LTE are by the line thick- 
ness (thick and thin lines, respectively). As a result, four profiles 
are shown for each of the 9 lines: iJ^ ILTE (thick solid line) _Rq' lte 
(thick dashed line), Rj; TE (thin solid line), and _R^ TE (thin dashed 
line). Shown in the ordinate is the normalized flux divided by the 
theoretical (pure) continuum; therefore, the local continuum level 
sometimes turns out to be slightly less than unity because of the 
extended wings of H lines. (Note that the theoretical equivalent 
widths corresponding to ij^ LTE profiles discussed in Ap- 

pendix B were calculated with respect to the local continuum, i.e. 
the maximum level in the neighborhood of the line profiles, irre- 
spective of the scales in the ordinate.) (a) Na I 5890, (b) Na I 
5896, (c) Na I 8183, (d) Na I 8195, (c) Na I 6154, (f) Na I 6161, 
(g) Na I 5683, (h) Na I 5688, (i) K I 7665, (j) K I 7699, and (k) 
Li I 6708. 



